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In mean field approximation, we study a ciiiral soliton of the linear sigma model with two flavors 
at finite temperature and density. The stable soliton solutions are calculated with some appropriate 
boundary conditions. Energy and radius of the soliton are determined in a hot medium of constituent 
quarks. It is found that for T < Tc, the energy of the soliton E* is less than the energy of three 
constituent quarks 3Mg, but with the increasing of temperature, the difference between E* and 
3Mq becomes smaller and smaller, once T > Tc, there is a sharp delocalization phase transition 
from hadron matter to quark matter coincident with the restoration of chiral symmetry. In the 
transition region, the thermodynamic properties show large discontinuities which is an indication 



, for a first-order phase transition. 

Ch ■ PACS numbers: 12.39.Fe,12.39.Ba,12.38.Aw,11.10.Wx 



I. INTRODUCTION 



^ . . . 

O |. It is widely believed that at sufficiently high temperatures and densities there is a quantum chromo dynamics 
JL ' (QCD) phase transition between normal nuclear matter and quark-gluon plasma (QGP), where quarks and gluons 
^ ', are no longer confined in hadrons ||] . Experimentally, the study of the QCD phase transition is supported by the 
(~| ' heavy-ion collisions in laboratories at ultrarelativistic energies, such as the Relativistic Heavy-Ion Collider (RHIC) 
at Brookhaven National Laboratory and the Large Hadron Collider (LHC) at CERN. These conducted experiments 
provide us with a chance to create a hot QCD matter and elucidate its properties. In order to explore of a wider range 
^ I of the QCD phase transition up to several times of the normal nuclear matter density, the new Facility for Antiproton 
. and Ion Research (FAIR) at Darmstadt, the Nuclotron-based Ion Collider Facility (NICA) at the Joint Institute for 
Ps) ' Nuclear Research (JINR) in Dubna and the Japan Proton Accelerator Research Complex (J-PARC) at JAEA and 
, KEK, will make such extreme conditions possible through collisions Q. 

On the theoretical side, due to the property of confinement, in nonperturbative regime of large distances, or 
equivalently low energies, the analytical as well as numerical methods have not been developed enough to fulfill the 
■ solution of low-energy nonperturbative cases, especially if the baryons are involved. Therefore the challenge to nuclear 
I physicists is to find models which can bridge the gap between the fundamental theory and our wealth of knowledge 
about low energy phenomenology, and these models should be successful in explaining empirical facts at low energies, 
for example, the dynamical breaking of chiral symmetry and the confinement, which arc both intimately related to 
1 J , the nonperturbative structure of the QCD vacuum. 
rN ' The linear sigma model (LSM) Q for the phenomenology of QCD has been proposed to describe the vacuum structure 
, with incorporating chiral symmetry and its spontaneous breaking. This model has two charming features according 
to two distinct low-energy phenomena. On the one hand, if we assume the thermal medium is mimicked by a uniform 
constituent quark medium with a dynamically generated mass, the model can be used to describe a restoration of 
the chiral symmetry at finite temperature and density in satisfactory agreement with the lattice calculations Q |^ and 
of the Nambu-Jona-Lasinio (NJL) model as well[0. On the other hand, starting from the same Lagrangian, bound 
states (chiral solitons) of valence quarks can be constructed through the interaction with a and tt mesons ^, and 
the nucleon naturally arises as a non-topological chiral soliton in the model. Moreover, the model has proven to 
be a successful approach to the description of nucleon static properties in vacuum]^, p^p^. Combining these two 
features together by requiring a soliton embedded in the hot soup of constituent quarks, the model seems to provide a 
suitable working scheme to simultaneously study both the restoration of chiral symmetry and the possible dissolution 
of the soliton, which simulates the delocalization (or deconfinement in some literatures) transition of nuclear matter 
to quark matter. Thus such a chiral soliton model is obviously more advantageous than other nontopological soliton 
modcls[|l5]|16 [1^ in the description of the hadron-quark phase transition||l^ jTof . 
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Using this model, the nucleon as a B = 1 chiral soliton in a cold quark medium has been investigated in Ref . In 
their studies, the nucleon is treated as a chiral soliton of the LSM, whereas the parameters /^r, and are chosen 
to be the medium-modified meson values within the NJL model. For finite temperature, more recently, by adopting 
the one-loop phenomenological mesonic potential and the coherent-pair approximation [|l3| , Abu-Shady and 
Mansour have studied nucleon properties at finite temperature Since the equally one- loop contribution of the 
quarks in potential is ignored, their studies can not be applied to finite density directly. In order to get a full 
picture of hadron-quark phase transition either at finite temperature or finite density, in present work, in mean field 
approximation, we are going to investigate the nucleon properties as well as the thermodynamics of the system at 
finite temperature and density. As we know the NJL model and the LSM treat the contribution of the Dirac sea 
differently, in the NJL model it is included explicitly up to a momentum cutoff A, while in the LSM this contribution 
is renormalized out. So that Different from the chiral soliton in the NJL model|2^]|Q, the present model has the 
benefit of renormalizability. Eventually, such a model would be investigated beyond mean field by the loop expansion, 
especially by using the CJT effective potential [p| ]25|. 

The structure of the paper is as follows: in the next section we introduce the chiral soliton model with two quark 
fiavors, nucleon static properties in the vacuum are briefly discussed and parameters are fixed. Then in mean field 
approximation, the chiral soliton solutions of the model at different temperatures and densities are self-consistently 
solved in Sect. III. Section IV is devoted to study static properties of nucleon at finite temperature and density, but 
we leave the study of the hadron-quark phase transition until section V. At the end, we give discussions and summary 
in Sec.VI. 



II. THE MODEL 



The chiral effective Lagrangian of the SU{2)b, x SU{2)l symmetry linear sigma model with two quark flavors has 
the formgH 

C = ^ [ijt^d^ - g{c7 + ij5T ■ Tf)] V + ^ {df^ad^cT + d^ir ■ d^'n) ~ C/(a, tt) , (1) 

here we have introduced a flavor-blind Yukawa coupling g of the isodoublet, spin-i quark flelds ip = {u, d) to interact 
with the spin-0, isosinglct a and the isotriplet pion field tt = (tti, 772, tts). The potential for the a and tt is parametrized 
as 

U{a, n) = ^ {a' + 7?^ - ^')' -H<J-'^ + ,Cml, (2) 
and the minimum energy occurs for chiral fields a and tt restricted to the chiral circle: 

'j'+^' = .fl (3) 

where = 93MeV is the pion decay constant and ~ 138 MeV is the pion mass, the last two constant terms in 
Eq.(H) are used to guarantee that the energy of vacuum in the absence of quarks is zero. The constant H is fixed by 
the PC AC relation which gives H = f 7^111^. Wc further assume that the chiral symmetry is spontaneously broken in 
the vacuum and the expectation values of the meson fields are (cr) = and (7?) = 0, then the dimensionless coupling 
constants g and A are the only two free parameters of the model, which are in turn conveniently re-expressed in terms 
of the constituent quark mass in vacuum, Mq = g/^, and the sigma mass, = + 2A/^. Finally, the quantity i?^ 
can subsequently be expressed as i?^ = — m\/\. 

The mass of the sigma meson is still a poorly known number, but the most recent compilation of the Particle Data 



Group considers that m„ can vary from 400 MeV to 550 MeV with full width 400 - 700 MeVl27j. The coupling 
constant g is usually fixed by the constituent quark mass in vacuum within the range of 300 ^ 500 MeV, which gives 
g ~ 3.3 ~ 5.3. In this work wc take rria = 472 MeV and g = 4.5 as the typical values in order to describe the 
properties of nucleon in vacuum successfully. 

In vacuum, the a and tt are taken as time- independent, classical c-number fields, which only differ from their 
vacuum values in the neighborhood of the quark sources. The state of the quarks {0„(r)} with energy {£«} and the 
(T(r), 7r(r) meson fields satisfy the coupled set of the Eulcr-Lagrange equations of motion 

- id ■ V(?!)„(r) - g/3[(j{r) + ij^f ■ n{r)](j)n{r) = £„(?!)„ (r), (4) 
_V2a(r) + = -g 0„(r)</>„(r) (5) 



do 

d-ff 



-V^7f{r) + = -<? E Mr)^l5fMr) (6) 
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with 

' 4'i{r)<Pn{v)d\^l (7) 



where d and /3 are the conventional Dirac matrices. 

The ground state of the chiral soliton is the state where N quarks in the same lowest Dirac state with energy 
e. In the following, our discussions are constrained in the case of = 3 for baryons. In order to obtain solutions of 
minimum energy, we adopt the "hedgehog" ansatz with the meson fields are spherically symmetric and valence quarks 
are in the lowest s-wave level 

cr — crir), tt = f7r(r), (8) 



u{r) 
ia ■ rv{r) 



X, (9) 



where x is a state in which the spin and isospin of the quark couple to zero: 

(a + f)x = 0. (10) 

Now the system is spherical symmetric and the Euler-Lagrange equations of motion (^-(^ transform in radial 
coordinates to 



duir) 
dr 



{e + ga{r))v{r) - gTr{r)u{r), (11) 



dvir) _ (2_ ^^^^^^ ^^^^ ^^^^^^ ^^^^^ ^^2^ 



dr \ r 



d^a{r) 2 dcr(r) dU 



dr"^ r dr da 
dPii{r) 2dTr{r) 27r(r) dU 



-^^Ng{u\r)~v\r)), (13) 
-2Ngu{r)v{r), (14) 



and the quark functions should satisfy the normalization condition 

47r J r^{u^{r) + v^{r))dr = 1. (15) 

These equations are subject to the boundary conditions which follow from the requirement of finite energy: 

w(0) =0,^^ = 0,^(0) =0 (16) 
dr 

m(oo) = O,cr(oo) :=:/^,7r(oo) = 0. (17) 

The asymptotic vacuum value of the soliton field has to be determined by an additional condition that the physical 
vacuum is recovered at infinity. In this "physical" vacuum the quarks are free Dirac particles of the constituent mass 
gai,, and the chiral symmetry is spontaneously broken. 

If we put N quarks into the lowest state with energy e, the total energy of the hedgehog baryon is given by 



= TVe + 47r 



1 f da\^ 1 f dn^ ^ 



2\dr 2 \dr 



Uia,7:) 



dr. (18) 



The model has two adjustable parameters g and A which can be chosen to fit various baryon properties, such as 
masses, charge radii and magnetic moments. Once the solutions to the above equations are obtained, one can calculate 
these physical quantities pertaining to the three-quark system, which have been measured experimentally. 



Unlike the Friedberg-Lee model[15 and its descendant models|17|, where the confinement of quarks is approximated 
through their interaction with the phcnomenological scalar field a which is introduced to describe the complicated 
nonperturbative features of the QCD vacuum. In many of these models, at large radius r, the a field assumes its 
vacuum value cr^, but at small r, the a field has a value close to the second minimum of the potential near zero. It 
means that in physical vacuum state, the quark mass is more than IGeV which makes it energetically unfavorable for 
the quark to exist freely, so that the effective heavy quarks have to be confined in hadron bags. Similar to the MIT 
bag model, they need to introduce the bag constant to make the baryon state stable. However, for the chiral soliton 
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FIG. 1: The total energy of system E and the energy of 3 free constituent quark SM, as a function of the Yukawa coupling 
g. For g < gc, 3Mq is lager than E, there is no stable bound state of quarks, however, for q > Qc, the soliton solution is stable, 
here gc = 4.24. 



in 





FIG. 2: The quark fields in relative unit and the a, n fields scaled with as function of the radius r for the parameters taken 
as = 472 MeV and g = 4.5. 



model, the interaction between the meson fields and the quarks is essential for the formation of a stable soliton, the 
state could be bound only when the total energy of system is lower than the energy of three free constituent quarks 
in system. Hence the coupling constant g or Mq is the significant parameter in the present model. In Fig.^ It is 
numerically shown that the critical value of coupling constant for mo- = 472 MeV is gc = 4.24, below that there is no 
stable bound state of quarks. 

Here after taking the set of parameters iria- = 472 MeV and g ~ 4.5, we plot the cr, tt and quark fields profiles in 
arbitrary unit as functions of radius r in Figj^ and calculate the properties of the nucleon, such as i? = 1170 MeV, 
i? = 0.877 fm, /ip = 0.319efm and gA/gv = 1-24. For comparing we adopt experimental values for the proton as 
follows|2|l: E = (Mp + Ma) 1085 MeV, Rp = 0.877 fm, fip = 0.294efm and gA/gv = 1-25. Therefore it is proven 
that this set of parameters can describe the properties of nucleon at zero temperature in a reasonable way. 
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III. MEAN FIELD APPROXIMATION 



A convenient framework of studying phase transitions is the thermal field theory. Within this framework, the finite 
temperature effective potential is an important and useful theoretical tool. In this section, in order to investigate the 
temperature and the chemical potential dependence of the chiral soliton, let us consider a spatially uniform system 
in thermodynamical equilibrium at temperature T and quark chemical potential fi. In general, the grand partition 
function reads 



Z = Trexp[-(?i-/i7\/')/T] 



(19) 



ly.T I ^ {ln[l + e-(^'-^)/^] + ln[l + e-(^^+^)/^]} . (21) 



where = i Jq^^ dt Jy (fix, V is the volume of the system and ^ = Mb/3. 

We evaluate the partition function in the mean- field approximation similar to the work of ||2^ . Thus we replace 
the meson fields by their expectation values in the action. In other words, we neglect both quantum and thermal 
fluctuations of the meson fields. The quarks and antiquarks are retained as quantum fields. The integration over the 
fermions yields a determinant which can be calculated by standard methods pSt. This generates an effective potential 
for the mesons. Finally, we obtain the thermodynamical potential density as 

f)(T,A.) = ^^ = f/(a,^) + r!^-;^, (20) 
with the quarks and antiquarks contribution 

Here, Uq ~ 2NfNc ~ 12 and Eg ~ ^Jjp + Af| is the valence quark and antiquark energy for u and d quarks, and 

the minus sign is the consequence of Fermi-Dirac statistics. On the chiral soliton background, the constituent quark 
(antiquark) mass Mq is defined as 

M'q=g^a'{r)+n'{r)). (22) 

As mentioned in the above discussion, as long as the coupling constant g is large enough, the minimum energy both 
in a vacuum and in an external hot/dense medium occurs for chiral fields restricted to the chiral circle: 

a^{r) + 7r'(r) = al (23) 

here the chiral order parameter is used to characterize the chiral symmetry breaking in space. In the following 
discussions, for a good approximation, the cr„ is taken as a constant for all r since they only slightly deviate from this 
constant within half of the radius R. Then the valence quark (antiquark) possesses almost same constituent mass 
in whole the space. This makes our discussions more reasonable than previous studies in the Friedberg-Lee model 
in thermal medium [p^ , where inside the soliton bag, the free quark carries the current mass, whereas outside the 
soliton bag, the confined quark has an unphysical large masses. Hereafter, a hot and dense thermal medium can be 
approximately modeled as a uniform constituent quark medium with solitons embedded in. 

The value of (7^ in equation ( p^ ) in thermal medium should be determined by minimizing the thermodynamical po- 
tential in Eq.(^) with respective to the classical meson fields cr(r) and 7r(r), therefore, in the mean-field approximation 
we arrive at a set of coupled equations for mesons 

d^a(r) 2 da(r) dfl , 9, . 9, 

—7V + -—r^-^=^9iu'{r)^v^{r)), (24) 
ar^ r ar oa 

d^Tr(r) 2dTT(r) 27r(r) dVL , ^ , ^ 

— ^ + -— ^ ^ - — = -27Vgu(r)z;(r), (25) 

where 

dVl 
9^ 



dVL 
d-K 



A (a2(r) + 7r2(r) - t?^) a{r) -H + gp,{r), (26) 
A {a^{r) + 7r2(r) - t?^) ^(r) + gpp,(r). (27) 
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The scalar and pseudoscalar densities of valence quarks and antiquarks can be expressed as 

f cPp 1 

Now the set of equations (|ll|) , , (^) and ( ^5|) can be solved for the soliton fields as functions of radius at finite 
temperature T and chemical potential fi under some appropriate boundary conditions, and the solutions of these 
coupled equations determine the behavior of the chiral order parameter cr„ as a function of T and fi. 

In the end of this section, it is worth pointing out that in our approach we have assumed a physical picture in which 
the valence quarks bearing a hedgehog configuration in the presence of the temperature and density, is simulated by 
these valence quarks facing the modified meson fields in thermal medium. So in calculation of the effective potential, 
the plane-wave valence quark states have been used, while, in fact, the quarks should be confined in the finite-size 
solitonic configuration. 



1 



1 



1 1 



(28) 
(29) 



IV. NUCLEON STATIC PROPERTIES AT FINITE TEMPERATURE AND DENSITY 

In order to gain the chiral soliton solutions and study nucleon static properties at finite temperature and density, 
the appropriate boundary conditions for the equations of motion in thermal medium need to be defined. As the case 
in vacuum, the a and tt are taken as time-independent, classical c-number fields, also we require they can only differ 
from their expectation values in the neighborhood of the quark sources. Therefore, in the absence of the quark source 
term for large r, the equations ( p^ and (|2^) reduce to the gap equations p9[| 

A {(J^r) + 7r2(r) - d^) a{r) ~H + gps{r) = 0, (30) 
A (a^ir) + 7r2(r) - 7r(r) + gpp,{r) = 0. (31) 

with spontaneously broken symmetry, <7{r) = (T„, and 7r(r) = 0. After self-consistently solving these gap equations 
( ^ ) and (^) with the constituent quark (antiquark) mass Aig in Eq.(]2^), we can get the value of the chiral order 
parameter (t„ for certain T and /i, to which the thermodynamical potential ^1{T, /i) has an absolute minimum. Now 
for the hedgehog baryon, we can safely define the boundary conditions for the coupled nonlinear equations of motion 
in above as 

v(0) = 0, = 0, 7r(0) = (32) 

ar 

u{oo) = 0, cr(c«) = (Jy, 7r(c«) = 0. (33) 

and these boundary conditions are indeed satisfied the requirement of finite energy of system, this in turn allow us to 
define the mass of nucleon (the total energy of the system) in thermal medium properly. 

In sequence, according to the boundary conditions Eqs.(^2|)(|3^), the set of equations of motion has been solved 
self-consistently by following the procedure as follows. For a given value of temperature T and density /i, the gap 
equations ( ^0[ ) and (^) can be numerically solved to obtain a constant chiral order parameter (jy which acts as the 
asymptotic value for a or a mass of the constituent quarks Mq. According to this particularly homogeneous solution 
(T„, the above coupled differential equations (pl|) , (|T^) , (p^ and ( p5| ) can be solved by the same method used in vacuum. 
After that we get the cr, tt, u and v field configurations as a function of distance radius r for certain T and p. In 
the end, by using these solutions, modifications of nucleon properties in thermal medium and the thermodynamics of 
system could be examined and investigated in detail. 

We firstly investigate soliton solutions at finite temperature and density. In Fig.^, we plot the u(r), w(r), u{r) and 
v{r) fields at zero and finite chemical potential {p ~ 200 MeV) for different temperatures. It is shown that all the 
fields are moving towards to the trivial zero values with temperature increasing. When T is lager than some critical 
temperature Tc, there only exist the trivial solutions for the coupled equations of motion and solitons are melted away. 
Furthermore, these trivial solutions indicate the restoration of the chiral symmetry in full space. The lack of solitonic 
solutions are taken as a signal for the delocalization of the baryonic phase. Sometimes, it is tempting to identify this 
with deconfinement; but such a conclusion seems to be out of the scope of the model itself, since the linear sigma 
model only incorporates the chiral symmetry, and according to the present model picture the nucleon is treated as 
bound state rather than the absolute confinement object. 
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FIG. 3: (a) The quark fields in relative unit and the a, n fields scaled with as function of the radius r at ^ = MeV, where 
the solid curves are for T = MeV, the dashed curves are for T = 100 MeV and the dash-dotted curves are for T = 116 MeV. 
(b) The quark fields in relative unit and the a, tt fields scaled with as function of the radius r a,t fj, = 200 MeV, where the 
solid curves are for T = MeV, the dashed curves are for T = 70 MeV and the dash-dotted curves are for T — 80 MeV. The 
parameters are taken as rricr ~ 472 MeV and g = 4.5 
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FIG. 4: The total energy of system E* and the energy of 3 free constituent quark 3Mq as function of the temperature T for 
the parameters taken as mcr ~ 472 MeV and g — 4.5. Here one set is for ^ — MeV, another set is for fi = 200 MeV and the 
meson interaction energy is not included both in vacuum and in thermal medium. 



Base on above analysis, as T is lower than Tc, there really exists the baryonic phase, but the stability of this phase 
should be checked carefully by comparing the total energy of the system in thermal medium with the energy of three 
free constituent quarks. By subtracting the homogeneous medium contribution [^3|, the total energy of system E* is 
given by the sum of the energy of the valence quarks and the kinetic energies of a and tt: 



E* = Ne + Att 



1 fda 

2 V 



1 fd-K 

2 V 



dr. 



(34) 



here, since we take the cr„ as a constant for all r, the meson interaction energy is discarded. This energy E* 
considered as an effective mass of nucleon in thermal medium. 
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FIG. 5: The proton charge r.m.s radius of a stable chiral sohton as a function of temperature T at /i = MeV and /i = 200 
MeV for the parameters taken as ttio- = 472 MeV and g = 4.5. The sohd curve is for = MeV while the dash-dotted curve 
for ^ = 200 MeV. 



In Fig.^, the total energy of system E* is plotted as a function of the temperature for = MeV and ^ ~ 200 
MeV, one can see that the total energy E* monotonically decreases with increasing temperature T from zero to 
higher value. As the temperature is close to the critical temperature Tc, E* starts to deviate from the ones in vacuum 
significantly, when T > Tc, E* jumps to zero quickly, which indicates the delocalization phase transition from nuclcon 
matter to quark matter. The energy of three constituent quark 3iVfq (or (t„) shows the similar behaviors as E* . The 
critical temperatures are 116 MeV and 80 MeV for /.t = MeV and ^ = 200 MeV, respectively. These results are 
very close to previous studies in the Friedberg-Lee model |l^ Q or the improved quark mass density-dependent 
model(IQMDD)|l^, where a first-order phase transition is predicted and the critical temperature at zero chemical 
potential is around 100 MeV. As discussed in Rcf. |2^, the critical temperature and the order of phase transition are 
strongly dependent on the coupling constant g. li g gets smaller, then Tc for the chiral phase transition will increase, on 
some critical value of g, the order of chiral phase transition would change from first-order to crossover. The qualitative 
behaviors of the soliton energy and the constituent quark mass changing with the temperature and density are also 

2^ , where they used the NJL model instead of the LSM to study chiral solitons[^ |Q. In 
contrast, our results are different from the work of Abu-Shady and Mansour, in that, by taking the different scheme, 
they have obtained the nucleon mass monotonically increased with temperature increasing and reached a very large 
value then slightly decreased, when they took the critical temperature around 161 MeV. This different result arises 
from the different treatment on the homogeneous thermal medium contribution, in present discussion, it is subtracted 
in order to study the stability of the single soliton by comparing with the energy of the three free constituent quarks. 
However, this contribution should be included if we want to investigate the overall thermodynamics properties of the 
soliton system, i.e. the internal energy, the entropy and free energy p^. 

Comparing the two energies in Fig.^, we can show that for T < Tc the nucleon bound sate is stable and 3Eq is 
larger than E*^ but the difference decreases with the increasing of temperature, and the two energies start to cross 
over at the critical temperature Tc, after Tc, because the chiral symmetry is restored and there is no soliton solution 
anymore. Hence it is concluded that a chiral phase transition together with a delocalization phase transition from 
hadron matter to quark matter is going to take place at same critical temperature. 

The proton charge r.m.s radii i? of a stable chiral soliton as a function of temperature for /.j = MeV and /.j = 200 
MeV are illustrated in Fig.^, it gives a signal of a swelling of the nucleon when temperature and density increase. In 
both cases, at low temperature the R increases slightly with the temperature increasing, as T approaching to Tc, R 
will sharply grow and disappear. Another interesting result displays in Fig.|^ is the maximal radius R at hight density 
gets reduced when compared with that at low density, this affect hints that the solitons are easily to overlap with 
each other at high density than the ones at low density when the T is close to the Tc. 



found in the Refs.p0]|23 
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FIG. 6: The sigma mass (solid line) and pion mass (dash-dotted line) as a function of temperature for fj, = MeV (right pair) 
and for /.t = 200 MeV (left pair). 



V. HADRON-QUARK PHASE TRANSITIONS 



In previous section we have obtained the effective nucleon mass at finite temperature and density, for convenience, we 
set Mjv = E* in the following discussions. In this section we will simply discuss how to calculate the thermodynamical 
variables of the system from hadron phase to quark phase in the chiral soliton model. 

The hadron and quark phase can be distinguished by empirical facts and phenomena at low energy. At low 
temperature and low baryon density, the hadronic phase exhibits a dynamical breaking of chiral symmetry and the 
confinement, and the baryon and meson act as the active degrees of freedom here. However, at very high temperature 
or baryon density, quarks and gluons will be set free to play the dominant roles in QGP. Based on the chiral soliton 
model, in hadron phase, the free quarks is not the ground state of strongly interacting matter, whereas three valence 
quarks will form the bound state of the nucleon, then in hadron phase we only have baryons and mesons. On the 
contrary, when T > Tc, the solitons are going to dissolve, then the hadronic phase will eventually evolve to quark 
phase. In the mean time, when incorporated with the Polyakov-loop, the model can be used to describe thermodynamic 
properties of QGP satisfactorily || HQ. 

Since there is no nucleon in quark phase, all information of the system is contained in the grand canonical potential 
which is given by fl in Eq.(pO[) with the restoration of chiral symmetry, and the pressure of system is directly given by 
Pqp = — r2„ji„(T, fi). However, in the hadronic phase, the situation becomes a little complicated because the system 
should be considered as a collection of nuclear matter (solitons) interacted through the self-consistent exchange of 
a, uj, p mesons. Moreover, in the spirit of the relativistic mean-field theory (RMF)|37|| or the quark meson coupling 



model (QMC)[|| 
nuclei extensively B 




I , one can investigate the medium modifications of nucleon properties in nuclear matter and finite 
Tl||^. This is out of the scope of our topic and left for future study. In following, we merely 
give an idealistical analysis on this topic by taking the hadronic phase as a noninteracting hadron gas composed of 
nucleons and tt, a mesons with the effective masses M^v, and M„ in thermal medium. Under this scenario, it is 
straightforward to write down the pressure of the system in terms of nucleons and mesons for the hadronic phase [B 



Php = 



^^T I ^|ln[l + e-(^"-^")/^]+ln[l + e-(^"+''~)/^]} 



(2^) 

.../Hlm-e-n}-../^ {inline--]}. 



(27r)3 



(35) 



Where vn 



4 for nucleon, 



- 3 for pion and i^o- = 1 for sigma meson. The energies Ej^ — \J "p^ ^ Mn'' 
^ are corresponding to nucleon, pion and sigma mesons, respectively. Here, Mi\ 



E^ ~ \/p'^ + AU'^, Ecr = Vp'^ + ^cr~ are corresponding to nucleon, pion and sigma mesons, respectively. Mere, Mat is 
obtained as the energy of soliton, whereas, the a and tt masses are determined by the curvature of ri(r, /i) in Eq.(|2C" 
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FIG. 7: The pressures as a function of temperature for different chemical potential. The solid curve is for fj, = MeV while 
the dash-dotted curve is for jj, = 200 MeV. 



at the global minimum: 



The sigma and pion masses for various T and ^ are shown in Figj6[ From this figure, the sigma mass first decreases 
smoothly and then rebounds and grows again at high temperature. The pion mass docs not change much at tem- 
peratures below Tc but then increases rapidly, approaching the sigma mass and signaling the restoration of chiral 
symmetry. At large temperature, the masses grow linearly with the temperature increasing. The sudden jumps in 
Fig.^ shows there is a first-order phase transition. 

Equipped with the masses for various temperature and densities, in Fig.^, we plot the pressures as a function of the 
temperature from hadron phase to quark phase at /x = MeV and ii = 200. From the figure, all curves show rapidly 
changed discontinuities at the critical temperature from hadron matter to quark matter. This indicates a first order 
phase transition for both the chiral phase transition and the delocalization transition, and signals a drastic structural 
changes in the system. Besides this, we can find that the Php in the hadronic phase is quite small when compared 
with the ones in quark phase, especially, for the case of higher baryon densities. This is due to the fact that the 
effective nucleon mass Afjv in our discussions slightly deviates from its vacuum value with temperature increasing. 
Only at the critical temperature Tc, experiences a sharp jump to zero. So that the contribution of the nucleon 
to the total pressure in Eq.(|35|) is too small in hadron phase. It subsequently leads to the calculated pressure in 
hadronic phase very different from that of the RMF theory at finite temperature and density[^|Q. Therefore, 
in order to confront the results in RMF theory in hadron phase, the model needs to be improved accordingly by 
including the interactions of nucleons and mesons with medium- modified coupling constants at hadron levelp^[p5[, 
or more practically, considering the possible overlap of nucleons. In the present work, however, we simply take this 
as a shortcoming of the simple model in the description of nuclear matter. 



VI. SUMMARY AND DISCUSSION 



We have studied the modification of the nucleon properties due to the restoration of the chiral symmetry at finite 
temperature and density within the linear sigma model with two flavor. The nucleon appears as a chiral soliton in the 
model, which is embedded in a thermal medium of constituent quarks with self-consistently determined effective mass. 
The chiral soliton solutions are solved in the mean-field approximation with the restriction to hedgehog configurations. 
The T and /i dependent energy of the single soliton is obtained. The only two free parameters g and rua are fixed in 
order to describe the properties of nucleon in vacuum successfully. The stabilities of the soliton solutions are analysed 
in thermal medium by comparing the effective mass of nucleon with the energy of three free constituent quarks. Our 
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results show that the chiral phase transition and the delocahzation phase transition from nucleon matter to quark 
matter take place simultaneously. For T < T^, the free constituent quarks are not the ground state of strongly 
interacting matter, the quarks will develop to form lower-energy bound states carrying the hedgehog configuration. 
However, as soon as the temperature T crosses over the Tc, such bound states can not live anymore, then the system 
experiences a first-order hadron-quark phase transition to the chirally symmetric phase. 

In this work, we predict Tc is about 116 MeV for = MeV, but 80 MeV for fj. = 200 MeV. It is much lower 
than the lattice predictions in the range of 150 — 200 McV. These relatively lower critical temperatures arc because 
of choosing a relative lager coupling constant g which has the functions of making the soliton stable as well as best 
fitting the proton charge r.m.s radius to 0.877 fm within the empirical parameter space of TOq-- To overcome this 
flaw, we can introduce the gluon interaction in the model at quark level, for example, the Polyakov-loop extension of 
the model[|4j. For roughly and qualitatively estimate, assuming the Polyakov loop variable $(r) and its Hermitian 
conjugate $(r) approach to their expectation values as t — >■ oo, then we can solve the coupled equations of motion 
again to get the updated behaviors of the as a function T and /x. By taking the Polyakov loop potential in Ref.Q 
with To = 190 MeV, the critical temperature at zero chemical potential now would be about 160 MeV, it is then quite 
rational to confront the constraint of the lattice data. 

It has turned out that the description of the hadron phase as a non-interacting hadron gas of the nuclcons and 
mesons with medium-modified masses has underestimated the important effects of their interactions at hadron level, 
and these interactions should be introduced to further reduce the effective masses of nucleons in nuclear matter 
significantly. Otherwise, such a simple model can not be used to describe the properties of nuclear matter and finite 
nuclei. Therefore, it is of interest to improve this model to make it suitable for the complete study of the hadron-quark 
phase transition in the whole region. All these works arc in progress. 
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